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There is conclusive evidence that the dominant component of the matter density in 
the universe is dark. The most striking indication of the existence of dark matter (DM) 
is the observations of flat rotation curves for spiral galaxies indicating that DM in 
galactic halos is about ten times more abundant than the luminous component. Dynam- 
ical evidence for DM in clusters of galaxies is also compelling. In terms of the critical 
density, pc = SH^M^JStc with Hq = 100/ikmsec"^Mpc"^ and Mpi the Planck mass, the 
amount of DM inferred from dynamics of clusters of galaxies is ^dm = Pdm/pc ~ 0.3. 
In addition, the most natural inflation models predict a flat universe, i.e., Qq = 1, while 
standard big-bang nucleosynthesis implies that ordinary baryonic matter can contribute 
at most 10% to Qq. This means that about 90% of the matter in our universe may be 
dark. 

It is usually assumed that DM consists of a species of a new, yet undiscovered, massive 
particle we denote as X. It is also often assumed that the DM is a thermal relic, i.e., it 
was in chemical equilibrium in the early universe. 

The simple assumption that the DM is a thermal relic is surprisingly restrictive. The 
limit Qx ~ 1 implies that the mass of a DM relic must be less than about 500 TeV 0. 
This upper bound turns out to be fatal to the proposal that DM consists of charged 
massive particles (CHAMPs — C^) The present experimental limits on superheavy 
hydrogen and C~p atoms are compatible with the CHAMP scenario only if they are more 
massive than about 10^ TeV |^. Similarly, current limits from underground detectors 
exclude the possibility that halo DM consists of colored particles of mass less than 500 
TeV. The standard lore is that the hunt for DM should concentrate on particles with 
mass of the order of the weak scale and with interaction with ordinary matter on the 
scale of the weak force. This has been the driving force behind the vast effort in DM 
detectors. 

In view of the unitarity argument, in order to consider thermal supermassive dark 
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matter ,0 one must invoke, for example, late-time entropy production to dilute the abun- 
dance of these supermassive particles 0], rendering the scenario unattractive. 

In this Letter we argue that recent developments in understanding how matter is 
created in the early universe suggests the possibility that that DM in the Universe might 
be naturally composed of nonthermal supermassive states. The supermassive dark matter 
(SDM) particles X may have a mass possibly as large as the Grand Unified Theory (GUT) 
scale. We suggest a number of cosmological production mechanisms for nonthermal 
supermassive dark matter. It is very intriguing that our considerations resurrect the 
possibility that the dark matter might be charged or even strongly interacting! 

We discuss four production mechanisms. We first propose production during reheat- 
ing of the universe after inflation. We point out that if the reheat temperature is denoted 
as Trh, the present abundance of SDM is proportional to (2000Mx/TR/f)~'', rather than 
exp(— Mx/Tr//) as one might naively expect. We then suggest the possibility of SDM 
production in preheating, making use of previous work that considered production of 
massive particles for baryogenesis. We then review the possibility of gravitational pro- 
duction of SDM at the end of the inflationary era. Finally, we propose that SDM might 
be created in the collisions of vacuum bubbles in a first-order phase transition. 

There are two necessary conditions for an SDM scenario. First, the SDM must be 

stable, or at least have a lifetime greater than the age of the universe. This may result 

from, for instance, supersymmetric theories where the breaking of supersymmetry is 

communicated to ordinary sparticles via the usual gauge forces 0. In particular, the 

secluded and the messenger sectors often have accidental symmetries analogous to baryon 

number. This means that the lightest particle in those sectors might be stable and very 

massive if supersymmetry is broken at a large scale [0]. Other natural candidates for 

supermassive DM arise in theories with discrete gauge symmetries Q and in string theory 
^In this paper, supermassive implies much more massive than the weak scale (about 100 GeV). 
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and M theory In the M-theory case, stable or metastable bound states of matter in 
the hidden sector, called cryptons, seem to be favoured over other possible candidates 
in string or M theory, such as the Kaluza-Klein states associated with extra dimensions. 
A specific string model that predicts cryptons as hidden-sector bound states weighing 
about 10^^ GeV is exhibited in P|. 

The second condition for SDM is that the particle must not have been in equilibrium 
when it froze out {i.e., it is not a thermal relic), otherwise fix would be larger than 
one.0 A sufficient condition for nonequilibrium is that the annihilation rate (per particle) 
must be smaller than the expansion rate: ncr|f | < H, where n is the number density, 
a\v\ is the annihilation rate times the M0ller flux factor, and H is the expansion rate. 
Conversely, if the SDM was created at some temperature and Qx < 1, then it is 
easy to show that it could not have attained equilibrium. To see this, assume X's were 
created in a radiation-dominated universe at temperature T*. Then Qx is given by Qx = 
Q.-^{T^ /TQ)mxnx{T^) / p-yiT^) , where Tq is the present temperature. (In this paper we will 
ignore dimensionless factors of order unity.) Using the fact that p-fiT^) = H{T^)MpiT^, 
we find nx{T^) / H{T^) = {Qx/fl'y)ToMpiT^,/Mx. Since we may safely take the limit 
a\v\ < M^^, nx{T^)(y\v\/H{%) must be less than {VLx/ft^)TQMpiTjMl. Thus, the 
requirement for nonequilibrium is 



This implies that if a nonrelativistic particle with Mx ~ 200 TeV was created at < Mx 

with a density low enough to result in Qx ~ 1, then its abundance must have been so 

small that it never attained equilibrium. Therefore, if there is some way to create SDM 

in the correct abundance to give Qx ~ 1, nonequilibrium is guaranteed. 

An attractive origin for SDM is during the defrosting phase after infiation. It is 
•^In this paper fli refers to the present value for species i 




(1) 
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important to realize that it is not necessary to convert a significant fraction of tlie 
available energy into massive particles; in fact, it must be an infinitesimal amount. If a 
fraction e of the available energy density is in the form of a massive, stable X particle, 
then VLx = ^^-yiTnu /Tq), where T^h is the "reheat" temperature. For Vtx = 1, this 
leads to the limit e ^ 10^^''(10^ GeV/TR/^). We will discuss how particles of mass much 
greater than Trh niay be created after inflation. 

In one extreme we might assume that the vacuum energy of inflation is immediately 
converted to radiation resulting in a reheat temperature Trh- In this case VLx can 
be calculated by integrating the Boltzmann equation with initial condition Nx = at 
T = Trh- One expects the X density to be suppressed by exp{—2Mx /Tjih)', indeed, one 
finds Qx ~ 1 for Mx/T^h ~ 25 + 0.5 ln(m|^((T|f|)), in agreement with previous estimates 
ig that for Trh ~ lO^GeV, the SDM mass would be about 2.5 x lO^^GeV. 

A second (and more plausible) scenario is that reheating is not instantaneous, but 
is the result of the decay of the infiaton field. In this approach the radiation is pro- 
duced as the infiaton decays. The SDM density is found by solving the coupled sys- 
tem of equations for the infiaton field energy, the radiation density, and the SDM 
mass density. The calculation has been recently reported in Ref. with the result 
VLx ~ i^'x{o'\v\){2QQQTjiH / MxY ■ Note that the suppression in Mx/Trh is not exponen- 
tial, but a power law (albeit a large power). Another crucial feature is the rather large 
factor of 10^. This implies that for a reheat temperature as low as lO^GeV, a particle of 
mass lO^^GeV can be produced in sufficient abundance to give Vtx ~ 1. 

The large difference in SDM masses in the two reheating scenarios arises because 
the peak temperature is much larger in the second scenario, even with identical Tjm. 
Because the temperature decreases as a~^/^ (a is the scale factor) during most of the 
reheating period in the second scenario, it must have once been much greater than Tr^. 
If we assume the radiation spectrum did not depart grossly from thermal, the effective 
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temperature having once been larger than Trh imphes that the density of particles 
with enough energy to create SDM was larger. Denoting as T2 the maximum effective 
temperature for the second scenario, T^jT^u ~ (M^/F^)^/"^ ^ 1, where V ^ is the effective 



decay rate of the inflaton. See [O for details. 

Another way to produce SDM after inflation is in a preliminary stage of reheating 
called "preheating" where nonperturbative quantum effects may lead to an ex- 



tremely effective dissipational dynamics and explosive particle production. Particles can 
be created in a broad parametric resonance with a fraction of the energy stored in the 
form of coherent inflaton oscillations at the end of inflation released after only a dozen 
oscillation periods. A crucial observation for our discussion is that particles with mass 
up to 10^^ GeV may be created during preheating [|l^, |14|, |15|, and that their distribution 
is nonthermal. If these particles are stable, they may be good candidates for SDM. 

To study how the creation of SDM occurs in preheating, let's take the simplest 
chaotic inflation potential: V[<\)) = M|0^/2 with ~ 10^^ GeV. We assume that 
the interaction term between the SDM and the inflaton field is g'^(j)'^\X\'^. Quantum 
fluctuations of the X field with momentum k during preheating approximately obey 
the Mathieu equation, Xj^ + [A{k) - 2qcos2z]Xk = 0, where q = g^(j)^/AMl, A{k) = 
[k"^ + Mx)/M^ + 2q (primes denotes differentiation with respect to z = M^t). Particle 
production occurs above the line A = 2q m an instability strip of width scaling as 
gi/2 fQj. loj-gQ g xhe condition for broad resonance, A — 2q < g^/^ [|T2| , p!3| , becomes 
(A;2 + M'j^)/Ml < g(t)/M^, which yields E\ = k"^ + Mj^ < gcpM^ for the typical energy 
of particles produced in preheating. Here is the amplitude of the oscillating inflaton 
field [O. The resulting estimate for the typical energy of particles at the end of the 



broad resonance regime for ~ lO-^Mpi is Ex ~ lO'^^^/yM^Mpi ~ g^l'^W^ GeV. 
Supermassive X bosons can be produced by the broad parametric resonance for Ex > 
Mxi which leads to the estimate that X production will be possible if Mx < g^^'^lQ^^ 



GeV. For (7^ ~ 1 one would have copious production of X particles as heavy as lO^^GeV, 
i.e., 100 times greater than the inflaton mass, which may be many orders of magnitude 
greater than the reheat temperature. In fact, in an expanding Universe and are 
time-dependent quantities and one should not only have very large field at the very 
beginning of the process of preheating, but also have E\ < g4>M^ until the end of 
preheating [16|. These considerations lead to an estimate of the upper bound on Mx 



slightly smaller than GeV |16|. Scatterings of X fluctuations off the zero mode 



of the inflaton field considerably limits the maximum magnitude of X fluctuations to be 
(^^)max ~ Mpg"^ 0. For example, (X^)^^^ < IQ-^^M^y if Mx = 10 M^. This restricts 
the corresponding number density of created X-particles. 

For a reheating temperature of the order of 100 GeV, the present abundance of SDM 
with mass Mx ~ 10^^ GeV is Qx ~ 1 if e ~ 10^^°. This small fraction corresponds to 
(X^) ~ lO^^^Mpi at the end of the preheating stage, a value naturally achieved for SDM 
masses in the GUT range [T7|. The creation of SDM through preheating and, therefore. 



the prediction of the present value of Qx, is very model dependent. The evolution of 
the background inflaton field responsible for the X production will be determined by its 
coupling to other fields since only a negligible fraction of its energy can go into SDM. 
We feel very encouraged, however, that it is possible to produce supermassive particles 
during preheating that are as massive as IO^^Tr//. 

Another possibility which has been recently investigated is the production of very 



massive particles by gravitational mechanisms [|T8], |19|. In particular, the desired abun- 
dance of SDM may be generated during the transition from the inflationary phase to 
a matter/radiation dominated phase as the result of the expansion of the background 



spacetime acting on vacuum fluctuations of the dark matter field [0. A crucial aspect 



of inflationary scenarios is the generation of density perturbations. A related effect, 
which does not seem to have attracted much attention, is the possibility of producing 
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matter fields due to the rapid change in the evolution of the scale factor around the 
end of inflation. Contrary to the first effect, the second one contributes to the homo- 
geneous background energy density that drives the cosmic expansion, and is essentially 
the familiar "particle production" effect of relativistic field theory in external fields. 

Very massive particles may be created in a nonthermal state in sufficient abundance 
for critical density today by classical gravitational effect on the vacuum state at the end 
of inflation. Mechanically, the particle creation scenario is similar to the inflationary 
generation of gravitational perturbations that seed the formation of large scale struc- 
tures. However, the quantum generation of energy density fluctuations from inflation is 
associated with the inflaton field which dominated the mass density of the universe, and 
not a generic, sub-dominant scalar field. 



If 0.04 < Mx/Hf. < 2 |18|, where Hf, is the Hubble constant at the end of inflation, 
DM produced gravitationally can have a density today of the order of the critical den- 
sity. This result is quite robust with respect to the fine details of the transition between 
the inflationary phase and the matter-dominated phase. The only requirement is that 
(iJe/10~®Mpi)^(TR///10^GeV) > 10^^. The observation of the cosmic background radi- 



ation anisotropy does not fix uniquely ife? but using Tjm < y/MpiH^, we find that the 
mechanism is effective only when > lO^GeV (or, Mx ~ lO^GeV). 



The distinguishing feature of this mechanism |118[ is the capability of generating 
particles with mass of the order of the inflaton mass even when the SDM only interacts 
extremely weakly (or not at all!) with other particles, including the inflaton. This feature 
makes the gravitational production mechanism quite model independent and, therefore, 
more appealing to us than the one occurring at preheating. 

Supermassive particles may also be produced in theories where inflation is completed 
by a first-order phase transition [^]. In these scenarios, the universe decays from its 



false vacuum state by bubble nucleation E^. When bubbles form, the energy of the 



false vacuum is entirely transformed into potential energy in the bubble walls, but as the 
bubbles expand, more and more of their energy becomes kinetic and the walls become 
highly relativistic. Eventually the bubble walls collide. 

During collisions, the walls oscillate through each other P2| and the kinetic energy 



is dispersed into low-energy scalar waves If these soft scalar quanta carry 



quantum numbers associated with some spontaneously broken symmetry, they may even 
lead to the phenomenon of nonthermal symmetry restoration [^]. We are, however, 
more interested in the fate of the potential energy of the walls, Mp = ArcqR'^, where rj is 
the energy per unit area of the bubble with radius R. The bubble walls can be imagined 
as a coherent state of inflaton particles, so that the typical energy E of the products of 
their decays is simply the inverse thickness of the wall, E ~ A^^. If the bubble walls 
are highly relativistic when they collide, there is the possibility of quantum production 
of nonthermal particles with mass well above the mass of the inflaton field, up to energy 
A^^ = 'jM^, 7 being the relativistic Lorentz factor. 

Suppose now that the SDM is some fermionic degree of freedom X and that it couples 
to the inflaton field by the Yukawa coupling g(f)XX. One can treat (the bubbles or 
walls) as a classical, external field and the SDM as a quantum field in the presence of this 
source. We are thus ignoring the backreaction of particle production on the evolution of 
the walls, but this is certainly a good approximation in our case. The number of SDM 
particles created in the collisions from the wall's potential energy is Nx ~ fxMp/Mx, 
where fx parametrizes the fraction of the primary decay products that are supermassive 
DM particles. The fraction fx will depend in general on the masses and the couplings of 
a particular theory in question. For the Yukawa coupling g, it is fx — S'^ln {'^M^/2Mx) 
p3| , p5[] . Supermassive particles in bubble collisions are produced out of equilibrium 
and they never attain chemical equilibrium. Assuming TpH — 100 GeV, the present 
abundance of SDM is Vtx ~ 1 if ~ lO^^a^/^. Here a^^ -C 1 denotes the fraction 
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of the bubble energy at nucleation which has remained in the form of potential energy 
at the time of collision. This simple analysis indicates that the correct magnitude for 
the abundance of X particles may be naturally obtained in the process of reheating in 
theories where inflation is terminated by bubble nucleation. 

In conclusion, we have suggested that DM may be elementary supermassive dark 
matter. Its mass may greatly exceed the electroweak scale, perhaps as the GUT scale. 
This is possible because the SDM was created in a nonthermal state and has never 
reached chemical equilibrium, thus avoiding the unitarity upper bound of about 10^ TeV. 
We have reviewed a number of ways SDM may be created. If reheating after inflation is 
preceded by a preheating stage, it is certainly possible to produce by resonance effects 
copious amounts of dark matter particles much heavier than the infiaton mass. We have 
also argued that the same may occur if inflation is completed by a first-order phase 
transition. These two scenarios are based on several assumptions about the structure of 
the theory, coupling constants, and the reheating temperature, but it is comfortable that 
the desirable abundance of nonthermal massive relics may be generated. Nonthermal 
SDM may be also created gravitationally at the end of inflation, with a significant mass 
range for which the SDM particles will have critical density today regardless of the fine 
details of the inflation- matter/radiation transition. This production mechanism involves 
the dynamics between the classical gravitational field and a quantum field; it needs no 
fine tuning of field couphngs or any coupling to the infiaton field. We are excited that 
the recent developments in understanding how matter is created in the early universe 
suggests that DM might be supermassive and resurrect the possibility that it might be 
charged or even strongly interacting. The implications of these fascinating options will 
be discussed elsewhere p6| . 
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